Adiabatic Landau-Zener-Stiickelberg transition with or without dissipation in low spin 

molecular system V15 
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The spin one lialf molecular system V15 shows no barrier against spin reversal. This makes 
possible direct phonon activation between the two levels. By tuning the field sweeping rate and the 
thermal coupling between sample and thermal reservoir we have control over the phonon-bottleneck 
phenomena previously reported in this system. We demonstrate adiabatic motion of molecule spins 
in time dependent magnetic fields and with different thermal couphng to the cryostat bath. We also 
discuss the origin of the zero-field tunneling splitting for a half-integer spin. 
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Research on quantum mechanical behavior in large sys- 
tems is of fundamental interest and also of major im- 
portance in quantum computation related applications 
([Q and refs. therein). A special class are the mag- 
netic two-level systems, realized in mesoscopic molecular 
crystals at miliKelvin temperatures Here, one scales 
up the quantum properties of one molecule in its local 
environment to large enough signals generated by the 
whole crystal. Any coupling between the spin-up and 
spin-down orientation of the molecular spin results in 
tunneling splittings whose magnitude control the quan- 
tum dynamics of the spin. Quantum tunneling of the 
magnetization across anisotropy barrier was observed in 
Mni2 or Feg high spin molecules j|] and quantum inter- 
ference of tunneling paths was demonstrated in Feg [||. 
In these molecules splittings are very small 10^^ K). 
Large splittings 10"'^ K) are usually found in low 
spin molecules or high and intermediate spin molecules 
in large transverse fields ^, ||. In such systems, the 
archetype of which is the low spin molecule V15 |^, ||] 
Landau-Zener-Stiickelberg (LZS) tunneling probability 
1^ is very close to 1 and, in the absence of environment, 
the transitions should be adiabatic. At finite tempera- 
ture, dynamics of the magnetic moment depend on direct 
spin-phonon transitions fed by the energy exchange with 
the cryostat. The S = 1/2, V15 molecule shows an in- 
trinsic phonon bottleneck (PB) effect with characteristic 
"butterfly" hysteresis loops. Such systems are experi- 
mental realizations of dissipative two- level systems . 
Similar results were reproduced since then in other low- 
spin systems |11 . The PB phenomena belong to a gen- 



eral class of almost adiabatic transitions in which only a 
small fraction of heat inflows during the process of ap- 
proach to equilibrium (which remains not completed). A 
similar example of almost adiabatic transition is found 
in the Foehn effect in meteorology and similar effects 
("magnetic Foehn" effect iQ) have been found in the 
magnetization processes of some Fe-rings . 

In this letter we present new results on mechanisms 



involved in the quantum dynamics of the V15 molecular 
spin. We prove fully adiabatical motion of this molecu- 
lar multispin system in time dependent magnetic fields 
and give more insight on the PB effect in showing ways 
to control dissipation in regions close to tunneling gaps. 
This gives the possibility to increase both relaxation and 
decoherence times. Furthermore, we discuss the origin of 
the zero-field splitting Aq for half-integer spin . 

The Vi5 complex (K6[Vi^As6042(H20)]8H20) forms 
a lattice with trigonal symmetry (a — 14.029 A, a = 
79.26°, V = 2632 A^) with two V15 molecules per cell 
||l||. AU the fifteen V^^ ions of spin 5 = 1/2 are placed 
in a quasi-spherical layered structure formed of a trian- 
gle sandwiched between two non-planar hexagons. The 
exchange couplings between the V ions are antiferro- 
magnetic making of V15 a typical example of frustrated 
molecule. After the PB study given in the next 
step, presented here, was to control the relaxation pro- 
cess by tuning the coupling between the spin-phonon sys- 
tem and the thermal reservoir. In order to reach a regime 
where spin dynamics is independent on the thermal bath, 
thus entering in a pure quantum regime, we minimized 
the coupling to the sample reservoir. This first limit of 
fast field variations and weak coupling to the cryostat, 
was reached by inserting a thermal isolator between the 
sample (a single-crystal of few tens of microns) and the 
sample-holder and restricting the temperature to lowest 
values. We optimized the other limit with best thermal 
coupling by using a mixture of silver powder with grease. 

The magnetization curves vs. applied field Bq (= 
fioH) of Fig. 1^ show either a weak hysteresis or no hys- 
teresis at all, depending on the competition between 
two field-dependent parameters: relaxation rate th and 
Aif/wAff, the time spent by system at a given two- 
level energy separation {Ah — ^Ag -t- {2^bBoY, vah = 
dAn/dt). In our experimental setup one may tune the 
ratio thvah/^h roughly between 10^"' — 10^ by adjust- 
ing the field sweeping rate and sample thermal coupling, 
which makes possible to study both regimes, spin sys- 
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FIG. 1: (a) Hysteresis loops for a sample rather well isolated 
from the cryostat (a — 130 sK^). The curve at 60 mK and 
0.14 T/s shows an adiabatic LZS spin transition, well fitted by 
Eq. ^. (b) Hysteresis loops measured in the case of a sample 
rather well coupled to the cryostat (a = 0.09 sK^). 



tem isolated (continous line, Fig. |^a) or in thermal equi- 
librium with the cryostat (dashed line. Fig. |^). The 
PB model accounts for the interpretation of such hys- 
teresis curves (for clarity only half cycles are shown in 
Figs. 1^, negative field half cycles being symmetric to the 
positive field ones). The spin temperature Ts is such 
that ni/n2 = exp(A///fcBTs') where ni_2 (fii,2eg) are 
levels out-of-equilibrium (equilibrium) occupation num- 
bers. Most molecules are out of equilibrium, at a tem- 
perature Ts different from the cryostat temperature T. 
Even for a good thermal contact between the sample and 
the cryostat, the energy flow from the latter is not suffi- 
cient to compensate the lack of phonons at energies hu = 
Ah ± Alj, if the field is swept fast enough. The phonon- 
bottleneck model presented in |l^ gives a relaxation 
law: -t/TH = x{t)-x{0) + \n{{x{t)-l)/{x{0)-l)) where 



TH 



tanh^(Aff/2/cBT) 



A2 



(1) 



is the relaxation time. The a parameter reflects the ther- 
mal coupling sample-cryostat (large value — weak cou- 
pling, small value — strong coupling); a is proportional 
to molecule density, phonon velocity, phonon-bath relax- 
ation time and inverse proportional to levels width (Aw, 
directly related to nuclear spin bath fluctuations ||l^ , of 
the order of 50 mK [0). To compare measured and cal- 
culated hysteresis curves one has to find the best set of 
parameters (a, Aq). The experimental curves in Fig. |l|a,b 
are best fitted by sets of parameters {a « 130 sK^, 
Ao « 80 mK) and (a « 0.09 sK^, Aq « 80 mK), re- 
spectively. As expected, {a « 0.15 sK^, Ao ^ 50 mK) 
obtained in shows a comparable Aq and a between 
the present ones, meaning that the coupling with cryo- 
stat was intermediate compared with those of Fig. |l|. It 
is interesting to mention that recent neutron scattering 
experiments done on V15 powder [jisl obtained a value of 
Aq comparable to our first determinations. 

In Fig. |l|b one can see that magnetization curve be- 
comes almost reversible if the field change is slow enough 
(dashed line). This is a natural consequence of Boltz- 
mann thermal equilibrium when spin-phonon transitions 
are highly probable. When sweeping rate vh or thermal 
isolation a are increasing, the PB phenomena becomes 
more present and manifests itself through an opening 
in the hysteresis loops. In very low fields (^ 0.1 T), 
the spin temperature is lower than the bath temperature 
{Ts < T). Then a PB plateau of almost constant mag- 
netization develops and Ts overpass T. At sufficiently 
high fields (~ 0.5 T) the system reaches its equilibrium 
{Ts = T) by a small phonon avalanche. Of a particular 
interest is the case a = 130 sK^ when the sample is ther- 
mally isolated from the sample holder. The reversibility 
of Fig. |l|a (continous line) is very different from the equi- 
librium reversibility. Here, the sweeping rate and isola- 
tion to the cryostat are large enough to drive the PB 
plateau near the saturation level: the PB phenomena 
is so important that the sample is completely isolated 
during our experimental time. The spin-phonon transi- 
tions have such a small probability that the system keeps 
its initial state (here, the groundstate) , allowing precise 
measurements of slow (spin-phonon transitions free) re- 
laxations. The multi-spin motion can therefore be con- 
sidered as adiabatic with an excellent approximation and 
magnetization curve should be given by: 



M = 



1 dA 



H 



2 dBn 




(2) 



Indeed, the adiabatic curve of Fig. |^ is nicely fitted by 
Eq. with Ao « 80 mK. 
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FIG. 2: Relaxation measurements at T=0.05 and 0.15 K 
(dots) fitted with the PB model (lines, see Eq. |l|), outside 
{B = 0.07 T) and inside the degeneracy zone (B = 0.014 T). 



Magnetic relaxation experiments were performed in 
the same conditions as in Fig. |l|a, for different values 
of the applied field. As an example, we show in Fig. ^ 
the results obtained for T ~ 0.05, 0.15 K and two fields 
especially chosen outside and inside the degeneracy zone 
(where Ah{Bq) is non- linear). The corresponding th val- 
ues, obtained by fitting the relaxation curves or simply 
calculated from Eq. |l| (with a = 130 sK^, Aq = 80 mK) 
are given in Tab. ^ . The fit using the relaxation law 
mentioned above (see Eq. |l|) is in remarkably good agree- 
ment with the data when the relaxation experiment is 
performed relatively far from the degeneracy point. The 
fit is less performant for the faster relaxations developing 
inside the degeneracy zone and the obtained value for th 
is 3 — 5 times shorter than the expected value (deduced 
from the general shape of the hysteresis curve). This 
shows that in this case, i.e. within the mixing region, the 
phonons bath is no longer sufficient to explain the dissipa- 
tion in this two level system. An excess relaxation comes 
from the Dzyaloshinsky-Moriya (D-M) interactions, al- 
lowed by molecule symmetry, and the spin bath, in par- 
ticular from the fast fluctuations of the ^^V nuclear spins. 
The D-M interactions generates off-diagonal terms cou- 
pling the spin-up and spin-down states, whereas the spin 
bath gives a noticeable spread of energy levels almost 
touching each other near zero-field, both effects enhanc- 
ing the relaxation in this particular region. Also, they 
are at the origin of the large zero-field splitting Aq in the 
Vi5 molecule as it will be discussed below. 

Vi5 has an effective spin 1/2 which is a half-integer. As 
well known as the Kramers theorem, all the energy levels 
have to be at least doubly degenerate in the time reversal 
symmetry. Therefore, if we consider the tunneling of the 
magnetization of the half-integer spins in a single mode 
path-integral formula, the tunneling rate must be zero. 
This is due to zero tunnel splitting Aq. The consequence 



TABLE I: Relaxation times obtained by fit of curves in Fig. C 
(3'''' column) and calculated (4*'' column) for a = 130 sK 
and Ao = 80 mK (obtained by fit of hysteresis loops Fig. Ha). 



is that magnetization reversals cannot be adiabatic with 
half-integer spins. However, as shown above, the V15 sys- 
tem shows very clear adiabatic spin reversal, with zero 
field splitting about almost 0.1 K. The most character- 
istic feature of V15 is its multispin character. We shall 
see that this character is a necessary but not sufficient 
condition for producing finite zero-field splittings with 
odd-integer spins. Note that Aq cannot be due to inter- 
nal dipolar fields which range in the mT scale. However, 
these fields together with nuclear spins certainly provide 
sources of decoherence. 

Due to frustration on antiferromagnetic triangles in 
Vi5, the energy structure of this system consists of one 
quartet {S = 3/2) and two doublets {S = 1/2) well sep- 
arated from excited levels by a gap of the order 10^ K. 
At zero field, the groundstate is fourfold degenerated (al- 
though with a spin S = 1/2). By making use of this 
degree of freedom one can have avoided level structures 
p9| . In order to mix the states whose magnetization dif- 
fers by 1, we need the interactions to contain and 
S~ . In the theoretical calculations the transverse field is 
often used for this purpose. However, in order to find the 
source of the avoided level crossing for the molecules in 
zero field, we have to find interactions with time rever- 
sal symmetry, that is, products of even number of spin 
operators, such as S^S'^, etc. Generally symmetric com- 
binations of operators, e.g. S^Sj + SfSj, are rearranged 
into anisotropy and do not cause avoided level structure. 
This is not the case with antisymmetric combinations of 
the form SfS^ - SfS^ for which it has been found that 
new avoided level structures are created (such terms were 
proposed to induce gaps otherwise forbidden in magnetic 
molecules [^). This is a sufficient condition for having 
a finite zero-field splitting with half-integer spins. The 
subsequent structure of energy spectrum at low energy 
is generally like the one depicted in Fig. ^ for the case 
of a frustrated triangle (the 15 spins 1/2 form 3 effec- 
tive spins 1/2, antiferromagnetically coupled |^). Here, 
we introduced the D-M interaction {D = 50 mK) in a 
simple way to demonstrate the gap actually opens with 
this type of perturbation. However, the real molecule 
consists of more complicated structure where a spin of 
the inner triangle connects to spins of upper and lower 
hexagons of the V15 molecule, and other D-M interac- 
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tions can exist and should be taken into account pi[ . 
More realistic energy structure will be investigated after 
the determination of arrangement of D-M interactions on 
the bonds [Q. For an effective triangle, there are two 
sets of avoided structures and they cross in zero field, al- 
lowing the Kramers theorem to be satisfied. Depending 
on the details of D-M interactions, two different situa- 
tions may occur depending if the two sets of level repul- 
sion structures are or are not orthogonal to each other. 
In the first case the four states are classified into two 
sets of subspaces, and the LZS mechanism exactly holds 
because the LZS transition occurs independently in each 
branch of avoided level crossing structure. In the second 
case, the four states form an irreductible space, where one 
of the state scatters into all of the four states when the 
field is swept, and the LZS mechanism should not hold. 
However, even in this last case, the change of magnetiza- 
tion at the crossing point as a function of the sweeping 
velocity is almost perfectly expressed by the LZS formula 
As a consequence, one can say that the change of 
magnetization in multispin systems with D-M interac- 
tions, should follows the LZS model in a wide range of 
interaction parameters. 

In short, previous studies of the "butterfly hysteresis 
loop" in the multi-spin V15 molecule Q are now ex- 
tended to the non-dissipative LZS regime. First, we have 
shown that the crossover from the dissipative to the non- 
dissipative LZS regime, can easily be achieved by chang- 
ing the coupling of the sample to the cryostat. This first 
demonstration of dissipation control in a two level sys- 
tem has important implications for further applications 
on quantum information and computers. The role of 
the spin-bath becomes visible when experiments are per- 
formed near zero-field, in the level mixing region. How- 
ever, the associated relaxation times are still long enough 
to allow adiabatic magnetization experiments. Secondly, 
we directly verified the zero- field gap previously obtained 
in the dissipative regime and attributed to the effect of 
D-M interactions of this multi-spin system. The third 
result is the clear comparison between experiments and 
theory, showing that the conditions for the existence of 
a gap in a half-integer multi-spin system with antisym- 
metric interactions {e.g. D-M interaction) are fuUfilled. 

We are very pleased to thank H. Bogge, A.K. Zvezdin, 
D. Mailly, K. Saito and N. Nagaosa for on-going collab- 
orations and useful discussions. 



* Present address: Department of Applied Physics, Delft 
University of Technology, 2628 CJ, Delft, Netherlands 

[1] A. J. Leggett, Suppl. Progress of Theoretical Physics 69, 
80 (1980), A. J. Leggett and A. Garg, Phys. Rev. Lett. 
54, 857 (1985), A. J. Leggett, Science 296, 861 (2002), 
M. N. Leuenberger and D. Loss, Nature 410, 789 (2001). 

[2] I. Tupitsyn and B. Barbara, in Magneto-Science- From 




-2 -10 1 

magnetic field (T) 



FIG. 3: Energy levels of 3 effective spins 1/2, obtained from 
H with J = 2.445 K yand D = 50 mK (for Aq = 80 mK). 



Molecules to Materials, edited by M. Drillon and J. Miller 
(Wiley VCH Verlag, 2001). 
[3] L. Thomas, F. Lionti, R. Ballon, D. Gatteschi, R. Sessoli, 
and B. Barbara, Nature 383, 145 (1996), J. R. Friedman, 
M. P. Sarachik, J. Tejada and R. Ziolo, Phys. Rev. Lett. 
76, 3830 (1996), C. Sangregorio, T. Ohm, C. Paulsen, 
R. SessoU and D. Gatteschi, Phys. Rev. Lett. 78, 4645 
(1997). 

[4] W. Wernsdorfer and R. Sessoli, Science 284, 133 (1999). 

[5] I. Chiorescu, R. Giraud, A. Jansen, A. Caneschi, and 
B. Barbara, Phys. Rev. Lett. 85, 4807 (2000), F. Luis, 
F.L. Mettes, J. Tejada, D. Gatteschi, L.J. de Jongh, 
Phys. Rev. Lett. 85, 4377 (2000), G. Bellessa, N. Vernier, 
B. Barbara, D. Gatteschi, Phys. Rev. Lett. 83, 416 
(1999). 

[6] I. Chiorescu, W. Wernsdorfer, A. Miiller, H. Bogge, and 

B. Barbara, Phys. Rev. Lett. 84, 3454 (2000). 
[7] I. Chiorescu, W. Wernsdorfer, B. Barbara, A. Miiller, and 

H. Bogge, J. Magn. Magn. Mat. 221(1-2), 103 (2000). 
[8] B. Barbara, I. Chiorescu, W. Wernsdorfer, H. Bogge, and 

A. Miiller, Progress Theor. Phys. Suppl. 145, in press 

(2002). 

[9] L. Landau, Phys. Z. Sowjetunion 2, 46 (1932), C. Zener, 
Proc. R. Soc. A 137, 696 (1932), E. C. G. Stiickelberg, 
Helv. Phys. Acta 5, 369 (1932), S. Miyashita, J. Phys. 
Soc. Jpn. 64, 3207 (1995), V.V. Dobrovitski, A.K. 
Zvezdin, EuroPhys. Lett. 38, 377 (1997). 

[10] A. Leggett et al. Rev. Mod. Phys. 59, 1 (1987), Y. Gefen, 
E. Ben- Jacob, A. O. Caldeira, Phys. Rev. B 36, 2770 
(1987), Y. Kayanuma, H. Nakayama, Phys. Rev. B 57, 
13099 (1998), M. Grifoni, P. Hanggi, Physics Reports 
304, 229 (1998). 

[11] R. Schenker, M. N. Leuenberger, G. Chaboussant, H. U. 
Giidel, and D. Loss, cond-mat/0204445. O. Waldmann, 
R. Koch, S. Schromm, P. Miiller, 1. Bernt, R. W. SaaL 
frank, cond-mat/0210309. 

[12] K. Saito and S. Miyashita, J. Phys. Soc. Jpn. 70, 3385 
(2001). 

[13] H. Nakano and S. Miyashita, J. Phys. Soc. Jpn. 70, 2151 
(2001), Y. Ajiro and Y. Inagaki, private communications; 



5 



Y. Narumi and K. Kindo, private communications. 
[14] A. Miiller and J. Doring, Angow. Chom., Intl. Ed. Engl. 

27(12), 1721 (1988), D. Gattcschi, L. Pardi, A.-L. Barra, 

A. Miiller and J. Doring, Nature 354, 465 (1991). 
[15] A. Abragam and B. Bleaney, Electronic paramagnetic 

resonance of transition ions (Clarendon Press, Oxford, 

1970). 

[16] N. Prokof'ev and P. Stamp, Rep. Prog. Phys. 63, 669 

(2000). 

[17] G. Carter, L. Bennett, and D. Kahan, in Metallic shifts in 
NMR, edited by B. Chalmers, J. Christian, and T. Mas- 



salski (Pergamon Press, London, 1977), vol. 20(1) of 

Progress in materials science. 
[18] C. Chaboussant ot al, cond-mat/0204365. 
[19] S. Miyashita and N. Nagaosa, Prog. Theor. Phys. 106, 

533 (2002). 

[20] B. Barbara, L. Thomas, F. Lionti, A. Sulpicc, and 
A. Caneschi, J. Magn. Magn. Mat. 177-181, 1324 (1998). 

[21] N. P. Konstantinidis and D. Coffey, cond-mat/0204435. 

[22] H. de Raedt, B. Barbara, I. Chiorescu, and S. Miyashita, 
in preparation. 



